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ABSTRACT 



Context. The design of modern instruments does not only imply thorough studies of instrumental effects but also a good understanding 
of the scientific analysis planned for the data. 

Aims. We investigate the reliability of Milne-Eddington (ME) inversions of high-resolution magnetograph measurements such as 
those to be obtained with the Imaging Magnetograph experiment (IMaX) aboard the Sunrise balloon. We also provide arguments to 
choose either Fe i 525.02 or 525.06 nm as the most suitable line for IMaX. 

Methods. We reproduce an IMaX observation using magnetoconvection simulations of the quiet Sun and synthesizing the four Stokes 
profiles emerging from them. The profiles are degraded by spatial and spectral resolution, noise, and limited wavelength sampling, 
just as real IMaX measurements. We invert these data and estimate the uncertainties in the retrieved physical parameters caused by 
the ME approximation and the spectral sampling. 

Results. It is possible to infer the magnetic field strength, inclination, azimuth, and line-of-sight velocity from standard IMaX 
measurements (4 Stokes parameters, 5 wavelength points, and a signal-to-noise ratio of 1000) applying ME inversions to any of 
the Fe i lines at 525 nm. We also find that telescope diff'raction has important eff'ects on the spectra coming from very high resolution 
observations of inhomogeneous atmospheres. Diffration reduces the amplitude of the polarization signals and changes the asymmetry 
of the Stokes profiles. 

Conclusions. The two Fei lines at 525 nm meet the scientific requirements of IMaX, but Fei 525.02 nm is to be preferred because it 
leads to smaller uncertainties in the retrieved parameters and off'ers a better detectability of the weakest (linear) polarization signals 
prevailing in the quiet Sun. 

Key words. Magnetic fields - Sun: photosphere - Instrumentation: high angular resolution 



1. Introduction 



With the advent of new instrumentation, we are witnessing a 
true revolution in our understanding of the quiet-Sun magnetism. 
The combination of superb spatial resolution and high-precision 
polarimetry has unveiled a panorama in which internetwork 
regions appear to be covered by very inclined, weak magnetic 
fields that had hitherto been elusive to observation. These 
findings have been possible thanks to an increase in the spatial 
resolution of spectropolarimetric measurements. Resolutions 
better than, say, 0'.'5 were hardly reachable with ground-based 
polarimeters just five years a go. Since 2006, however, the 
Hinode spectropolarimeter (SP: iLites et al.l l200ll: iKosugi et"aD 
I2007t It suneta et aP [20081), is providing data with spatial 
resolutions of 0'.'32 on a regular basis, with a signal-to-noise 
ratio (SNR) of about 1000. The analysis of these measurements 
has shown that the weak internetwor k fields occupy fractions 

[20081: 
Thus, 
form the 



instruments and their properties have to be determined by means 
of careful analyses. 

Further improvements in the spatial resolution will conceiv- 
ably lead to new discoveries. Since the polarization signals scale 
linearly with the filling factor, they will become larger if the in- 
ternetwork is unevenly permeated by magnetic structures not 
much smaller than the current resolution limit. On the con- 
trary, if the characteristic sizes are very small and the struc- 
tures homogeneously distributed, then the observed polariza- 
tion amplitudes will not increase at higher spatial resolution 
because the amount of cancellation of opposite polarity sig- 
nals will remain the sam43. There might be a third possibility 
that the magnetic structures are small but fractally distributed 
(iPietarila Graham et al. 1 1200^ . In this case, the cancellation of 
opposite polarity signals has been observed to depend on the spa- 
tial resolution at least down to the 0'.'3 of Hinod^ The subject 



of ~ 20 - 45% of the SP pi xel (iLites et al 
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200 



that 



iQrozco Suarez et al.ll2007l I201QI: [Orozco Suarez 
we are nearly resolving the magnetic fields 
internetwork, in the sense that the filling factor is approaching 
100%. But these fields are not yet directly visible with current 



^ Traditionally, the cancellation of Zeeman signals has been 
associated with the presence of mixed-polarity fields at the spatial 

ved by 

quiet Sun (see, e.g., Trujillo Bueno et al.i 2( 

^ This efi'ect was pre viously shown in simulations by 
ISanchez Almeida et all ( l2003h . 



]|2004). 
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Table 1. Basic parameters of the SuNRiSE/IMaX system. 

Telescope aperture 1 m 

Telescope central obscuration 35.2% 

Working wavelength 525.02 nm 

Spatial resolution ~ (X' 11 

CCD PIXEL SIZE (y!055xQf!055 

IMaX passband (FWHM) ~ 6 pm 



is currently feeding an intense debate. A distinction between the 
scenarios and the determination of their relative contributions to 
the total magnetic flux is very important and can efl'ectively be 
achieved by using larger telescopes and advanced techniques for 
high resolution spectropolarimetry. 

In moments of very good seeing, the CRisp Imaging Spectro- 
Polarimeter (CRISP; iScharmed 1200 6). installed at the Swedish 
Solar Telescope on La Palma (Spain), deliver polarimetric 
images near the difl'raction limit of its 1-m telescope with 
the help of adaptive optics and computatio nally-expensive 
image restoration techniques (Ivan Noort et al.l 2005). CRISP 
is a filter-based instrument, which allows the observation of 
large fields of view while preserving high spatial resolution 
and polarimetric accuracy. Compared to spectrograph-based 
polarimeters, filter instruments also have drawbacks like poorer 
spectral resolution and wavelength sampling. Several lines 
can be observed sequentially, but at the expense of slower 
cadences (without degrading the SNR). These advantages and 
disadventages imply diff'erences in the results obtained with 
each instrument that deserve attention, especially when a new 
one is designed and constructed. The way the instrument treats 
the polarization signals, the pointing stability, the wavelength 
samples, the magnetic sensitivity of the selected spectral line, the 
polarization modulation scheme, the polarimetric eflftciency, and 
the influence of environmental (instrumental plus atmospheric) 
polarization are some of the issues to be considered. 

The Imaging Magnetograph experiment (IMaX) is a dual- 
beam vector spectropolarimeter based on a Fabry-Perot etalon 
as the spectrum analyzer and a pair of li quid cry stal variable 
retarders as the pol arization mod i ilator (iMartinez P illet et al. 
i2010: see also J ochum e t al.ll2003l iMartinez Pillet et al. 2004 
and lAlvarez-Her rero et al.l l2006h IMaX is one of the post- 
focus instruments of the Sunrise Im telescope, which flew over 
the Artie in a stratospheric balloon from June 8 to 13, 2010 
(iBarthol et al.l l2010h . In this way, continuous observations of 
the Sun were possible avoiding most of the Earth atmosphere. 
Table [T] summarizes the basic parameters of the instrument. 
IMaX provides capabilities for obtaining polarimetric images 
near the diflraction limit of Sunrise (O'.'ll at 525 nm), with a 
SNR of about 1000. The instrument has several observing modes 
in which the number of wavelengths is varied from 3 through 
12 (one of them is always reserved for the continuum), but the 
"regular", so-called vector spectropolarimetric, mode is one that 
scans the line at 5 wavelength points in 33 s. Will it be possible 
to accurately infer the full vector magnetic field from the regular 
IMaX measurements using appropriate inversion techniques? 
Answering this question is the main aim of the present paper. 

To assess the accuracy to which the magnetic field strength, 
inclination, azimuth, and line-of-sight (LOS) velocity can be 
determined, we here simulate an IMaX observation. We use 
radiative magnetohydrodynamic (MHD) models to synthesize 
the Stokes profiles of the photospheric Fei lines at 525.02 and 
525.06 nm. These two lines, very close in wavelength and hence 
observable with the same pre-filter, were considered for the 



IMaX instrument since the beginning of the design phase. In 
fact, part of the results of this paper provided the rationale 
for selecting one of the lines instead of the other. An earlier 
investigation on magnetic inferences from filtergra ms with 
limited wavelength sampling was carried out by Graham et al.l 
(l2002h . They used Milne-Eddington (ME) profiles to simulate 
the observations. Our paper can be considered as an extension 
of their work, since MHD simulations provide a more realistic 
description of the solar photosphere. 

The simulated Stokes profiles are degraded by telescope 
diff'raction and detector pixel size to a spatial resolution of 
80 km on the solar surface, and then sampled at five wavelengths 
(four within the line plus another in the nearby continuum) to 
mimic the nominal vector spectropolarimetric mode of IMaX. 
From these data we determine the magnetic field vector and 
LOS velocity using ME inversions and compare them with the 
parameters resulting from the inversion of fully resolved and 
critically sampled profiles. To simulate a real observation, the 
typical noise of IMaX measurements is added to the profiles. 

The paper is structured as follows: Sect. [2] describes the 
spatial degradation of the simulated Stokes profiles and the 
eff'ects of the tunable filter on the spectra. In that Section we also 
analyze the eff'ects of noise and the drawbacks of selecting only 
a few wavelength positions across the line. Section [3] deals with 
the results of ME inversions. There we estimate the uncertainties 
of the retrieved physical parameters. In Sect. |4] we summarize the 
main conclusions of this work. 



2. Simulating IMaX observations 

2.1. Spatial degradation 

To simulate an IMaX observation we use model atmosphere s 
from the radiative MHD calculations of IVogler et al.l (l2005l) . 
They were initialized with a distribution of mixed polarity fields 
with an average magnetic field strength of (B) ~ 200 G at 
log T = -I. The snapshot used here was taken when (B) had 
decayed to about -140 G. We have chosen this particular run 
because it contains both network- and internetwork-like regions. 
The extent of the computational box is 288x288xlOO pixels, 
with a spatial sampling of 20.8 km (0'.'0287) in the horizontal 
direction and 14 km in the vertical direction. The bottom of the 
box is located 800 km below r = 1. Mo re details about this run 
can be found in lKhomenko et al.l (l2005l). 

Since the work of IVogler et al. I (l2005h . simulations of higher 
numerical resolution have become available. They certainly 
represent a significant step forward, but we believe they would 
not change our conclusions, at least for the larger magnetic 
structures present in the network and most of the internetwork 
- the ones observed by IMaX during its first flight and also 
by Hinode sin ce its launch in 20 06. The reason is that the 
simulations of IVogler et al.l (l2005b have been confronted with 
observations thoroughly with excellent results. This indicates 
that they already provide a very good description of the solar 
atmosphere. 

The MHD models are used to generate the observations 
by synthesizing the Stokes /, Q, U, and V profiles with the 
SIR code (Ruiz Cobo & D el Tbro Iniestal 119921) . The spectral 
synthesis is performed in a wavelength range that spans 1 nm 
and includes the Fei lines at 525.02 and 525.06 nm. The 
wavelength step is 1 pm. Next, the monochromatic images 
are spatially degraded as the SuNRiSE/IMaX system does. The 
modulation transfer function (MTF) of the whole system is 
depicted in Fig.[T]and includes the 1 m aperture of the telescope. 
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Fig. 1. Continuum intensity map at the original resolution of the MHD models (left) and degraded considering telescope diffraction 
and detector pixel size (right). The central panel shows the MTF of the CCD (dotted line), that of the diffraction-limited telescope 
including the central obscuration (dashed line), and the total MTF (solid line). 



the central obscuration of the entrance pupil (caused by the 
secondary mirror), and the effects of image pixelation. The 
spatially degraded images are binned by a factor 2, hence their 
sampling is almost identical to that of IMaX with a pixel size 
of (y.'055. The original continuum contrast of 17.9% decreases 
only by 4% due to diffraction and CCD pixelation. After that we 
convolve the Stokes profiles with a Gaussian of 6 pm FWHM 
to account for the limited spectral resolving power of the Fabry- 
Perot etalon. Finally, we add noise and select four wavelength 
samples across the line plus a wavelength point in the nearby 
continuum. 

The spatial degradation described above represents a rather 
pessimistic estimate because IMaX took in-flight calibrations 
of the optical aberrations for post-facto restorations. This will 
increase the power at all frequencies up to the cut- off limit. 
The calibration procedure (IVargas Dommguezl l2Q08h is based 
on the phase divers ity (PD) techni que for image restoration and 
wavefront sensing (iPaxman et al.l ll996). The instrument has a 
glass plate that can be optionally intercalated between the beam 
splitter and one of the CCDs to induce a controlled defocus. 
Therefore, a simultaneous focus-defocus image pair (PD pair) 
can be recorded by employing the two cameras. From this 
pair, an estimate of the wavefront aberration, parameterized in 
the form of a Zernike polynomials expansion, will be possible 
during the data processing by means of a PD inversion code. 
Deconvolution of the derived point spread function of the system 
will be carried out by Wiener filtering the data. Since this 
will partially correct the effects of the system, our using the 
theoretical MTF shown in Fig. 1 means we are considering a 
worst-case scenario. 

The two most important effects of telescope diffraction 
are spatial blurring and Stokes profile smearing. The first 
takes place because (polarized) light spreads from every 
pixel to its immediate surroundings. As a consequence, the 
spatial resolution decreases, the polarization amplitudes tend 
to diminish, and the signals are more vulnerable to photon 
noise. In addition, the spatial resolution decreases. Figure [2l 
displays monochromatic Stokes Q, U and V images at +7.7 pm 
from the center of the Fei 525.02 nm line where these effects 
can be observed. To emphasize details we represent a small 
area of the simulation snapshot. (A) stands for the images at 
original resolution and (B) for the images spatially degraded 
by telescope diffraction and CCD pixelation. We have added 



noise at the level of 10"^/c in order to see how it affects the 
observations. Since Stokes Q and U have weaker signals than 
Stokes y, they are significantly more degraded than the latter. In 
general, the images do not look strongly altered; the polarization 
amplitudes decrease and some structures escape detection, but 
the morphological aspect is essentially preserved. 

The second effect comes from the fact that the spatially 
averaged Stokes profiles are less asymmetric. As an example, we 
show in Fig. [3] the effect of pure spatial degradation (telescope 
plus detector) on the polarization profiles. The Stokes / and V 
signals from two different points (indicated by numbers 1 and 
2 in Fig. [T]) are displayed. Note the strong modification of the 
continuum and the shape of the profiles, prior to any action of 
the IMaX etalon. Similar modifications can be expected from 
any instrument, but apparently they have not been described in 
the literature. MHD simulations like the ones we are using here 
are the best means to date to bring this effect to the attention 
of the community: details of the profiles are lost because of the 
finite spatial resolving power. 

2.2. Spectral degradation 

IMaX uses a single Fabry-Perot etalon in double pass for the 
spectral analysis. The FWHM of the IMaX spectral point spread 
function (PSF) is ~6 pm. For simplicity, we approximate the 
shape of the PSF by a Gaussian function. The effect of the IMaX 
etalon on the polarization images is shown in column (C) of Fig. 
[21 Naturally, the amplitudes diminish and the weakest signals 
disappear below the noise. 

As is well known, the spectral smearing induced by the 
etalon smooths out the profile shapes and the Stokes spectra 
become more symmetric. Therefore, some of the weakest 
signals, typically linear polarization signals, cannot be detected 
after the combined action of the telescope, the etalon, and 
the detector. But the problem is not only on detection. These 
effects also reduce the information content of the detected 
profiles. The sensitivity to gas velocities, for example, decreases 
because t he etalon broadens the p rofiles and makes them 
shallower (Cabrera Sol ana et ani20Q5l) . The sensitivity to other 
atmospheric parameters also decreases with a wider spectral 
PSF because respon se functions also get broadened by it 
(iQrozco Su^z & Del Toro Iniestal[2QQ7l) . 
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Fig. 2. Monochromatic Stokes Q (top), U (middle), and V (bottom) images at +7.7 pm from the central wavelength of the Fei line 
at 525.02 nm. From left to right: the original data (A), the spatially degraded data (B), and the spatially and spectrally degraded data 
after taking into account the spectral PSF of 6 pm (C). Only a region of 150 x 150 pixels is displayed. 



2.3. The effect of noise 

Both telescope diffraction (plus CCD pixelation) and spectral 
analysis reduce the amplitude of the polarization signals. 
Consequently, they are more affected by the noise. In this section 
we quantify the influence of photon noise on the measurements. 

Figure|4]shows the percentage of pixels whose Stokes QorU 
amplitudes exceed specified noise levels. The Fe i lines at 525.06 
and 525.02 nm are represented in red and black, respectively. 
Obviously, the amount of detectable linear polarization signal 
decreases quite rapidly with noise. Noise levels of 3 x 10"^/c and 
10"^/c, the latter corresponding to the IMaX case, are marked 
with vertical lines. The dotted lines stand for the spatially 
degraded Stokes profiles and the solid lines for the spatially and 
spectrally degraded profies. 80% and 55% of the pixels show 
linear polarization signals in Fe i 525.02 nm and 525.06 nm that 
exceed the noise level before the light enters the instrument. 
These values change to 50% and 20% after the spectral analysis. 
Interestingly, the final detectability of linear polarization in the 
525.02 nm line is almost the same as that of 525.06 nm before the 
spectral analysis. Thus, Fe 1 525.02 nm presents clear advantages 
for detecting weak linear polarization signals in the quiet Sun. 



Figure[5] shows histograms of the Stokes V amplitudes in the 
original maps (solid lines), the spatially degraded maps (dashed 
lines), and the spatially plus spectrally degraded maps (dash- 
dotted lines), for Fei 525.06 and 525.02 nm (top and bottom, 
respectively). The histograms show an asymmetric distribution 
of amplitudes with two peaks corresponding to network and 
internetwork points, the latter being responsible for the smaller 
y amplitudes. The histograms for the 525.02 nm line are shifted 
to stronger signals because of the enhanced sensitivity of this 
line to magnetic fields. Again, specific noise levels are indicated 
by the vertical lines. In the original data, the percentage of pixels 
with Stokes V amplitudes below 10"^ h is 0.06% and 0.5% for 
the 525.02 and 525.06 nm lines, respectively. These percentages 
change to 1.6% and to 4.3% when the profiles are degraded 
spatially and spectrall}0. Therefore, most of the V signals would 
be observable with our instrument if we use the 525.02 nm 
line and the SNR is 1000. Thus, Fei 525.02 nm reveals itself 



^ Of course, if SNR = 300, the Stokes V detectability decreases: the 
percentages reach 2% and 5.3% in the original data, and 11.4% and 
23.2% after degradation. 
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Fig. 3. Stokes / and V profiles of Fei 525.02 nm corresponding 
to pixels #1 (top) and 2 (bottom) in Fig.[T] The solid and dashed 
lines represent the original simulations and the maps degraded 
by telescope diff'raction and CCD pixelation, respectively. The 
vertical lines indicate wavelength samples at [+8, ±4] pm from 
the central wavelength of the line. To make a pixel-to-pixel 
comparison, we have taken the Stokes profiles from the degraded 
data before binning the images to match the IMaX CCD pixel 
size. Therefore, dissimilarities are enhanced. 
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Fig. 5. Distribution of the Stokes V/Iq amplitudes in the 
original, spatially and spectrally degraded data. The vertical 
lines represent noise levels of 10"^ and 3 x 10"^/c. 

files were going to be obtained with the slit spectropolarimeter 
originally included in the scientific payload of Sunrise. Here we 
are concerned with the regular IMaX spectropolarimetric mode 
and therefore we have to acknowledge a further decrease in the 
fraction of polarization signals detectable above the noise level. 



0,000 0,001 0,002 0,005 0.004 0.005 
Noise (in units) 

Fig. 4. Variation in the percentage of pixels where either Stokes 
Q or U exceed a given noise level. The dotted and solid lines 
correspond to data before and after convolution with the spectral 
PSF. The vertical lines represent noise levels of 10"^ and 3 x 

as the better choice for detecting both weak linear and circular 
polarization signals. 

Further degradation is expected due to the limited wave- 
length sampling. As mentioned in the introduction, the regular 
vector spectropolarimetric mode of IMaX measures the polar- 
ization state in 4 wavelength within the line plus another in the 
continuum. These wavelengths have been marked with vertical 
lines in Fig. [3j The choice of a limited wavelength sampling is 
not arbitrary: on the one hand, scanning the line at just a few 
points reduces the scan time, improving the temporal integrity 
of the data and their immunity to variations in the conditions 
of the dynamic Sun or changes in the solar scene; on the other 
hand, a better sampling was not needed because full line pro- 



3. Inversion of the Stokes profiles 

In this section we address two questions. First, we want to 
know which line is less aff'ected by the constraints imposed 
by the instrument (i.e., spectral degradation, limited wavelength 
sampling, photon noise). The answer to this question is 
important to choose the most suitable line for IMaX. Second, we 
want to determine the uncertainties that can be expected from the 
analysis of the measurements taken in the selected line. 

To answer the first question, we need to evaluate the 
influence of instrumental efl'ects on the ME parameters resulting 
from the inversion of the two spectral lines. This requires 
a differential analysis. In principle, since the original model 
atmospheres are at our disposal, it would be straightforward to 
compare the inversion results for each pixel and each spectral 
line as observed by IMaX with the corresponding MHD models. 
However, such an absolute analysis would not inform us about 
the relative merits of the two lines. The reason is that they 
are formed in slightly difl'erent layers of the atmosphere, and a 
direct comparison of the ME parameters derived from them with 
the MHD parameters at a fixed height would yield diff'erences 
that do not reflect the capabilities of the lines, but rather their 
difl'erent formation heights. 
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The best way to solve this problem is to set up a reference for 
each of the lines. Here, the reference is taken to be the result of a 
ME inversion of the spatially degraded Stokes profiles without 
noise, sampled at 61 wavelength points separated by 1 pm. 
They represent the idealized case in which the instrument would 
measure fully resolved Stokes profiles with excellent (infinite) 
signal-to-noise ratios. Note that the reference includes the eff'ects 
of telescope diff'raction, whose influence on the polarization 
spectra has already been analyzed in Sect. O 

The differential analysis is performed by inverting the 
simulated IMaX measurements in the two lines and comparing 
the inferred ME parameters with the corresponding reference. 
The results, presented in Section [TT] led to the selection of Fe i 
525.50 nm as the IMaX line. 

Answering the second question requires a comparison of the 
original MHD models and the ME parameters derived from the 
Fe I 520.52 nm measurements. This is done in Sect. l3.2l below. 

All the inversions are carried out assuming a simple one- 
component model atmosphere and no stray-light contamination. 
The magnetic field vector is initialized with B = 200 G, y = 20°, 
and y - 20°. The inversion code used in the tests is MILOS 
(iQrozco Suarez & Del Toro Iniestall2007l) . 

3.1. Impact of instrumental effects 

To determine the deviations induced by the instrument with 
respect to the reference, we invert the simulated IMaX data, i.e., 
the spatially degraded Stokes profiles, smeared by a spectral PSF 
of 6 pm FWHM, sampled at 5 wavelength points (one of which 
is continuum), and with noise added at the level of 10"^ /c. 

The inversion returns the three components of the vector 
magnetic field and the LOS velocity. The maps corresponding 
to the Fei 525.02 nm line are displayed in Fig. [6j The left 
column shows the reference. The middle column displays the 
results of inverting the profiles convolved with a PSF of 6 pm 
FWHM and with noise, but still sampled at 61 wavelengths. The 
right column gives the results obtained from the simulated IMaX 
measurements. Only a region corresponding to 150 x 150 pixels 
is displayed. 

A comparison of the maps demonstrates that the field 
strength and the LOS velocity are well recovered. Only small 
deviations from the reference parameters can be seen, although 
the deviations increase as the field weakens (especially below 
100 G). The inclination and azimuth maps show the largest 
deviations because the spectral smearing, the noise, and the 
limited sampling reduce the Stokes Q and U amplitudes. The 
noisier areas in the maps derived from the simulated IMaX data 
correspond to regions where the linear polarization signals are 
buried in the noise and the field is very weak. 

The overall conclusion drawn from these tests is that ME 
inversions of IMaX measurements in the Fe i 525.02 nm line are 
reliable. The accuracy can be expected to be higher for the field 
strength and the LOS velocity, but also the magnetic inclination 
and azimuth, although noisier, will be recovered reasonably 
well. 

To support this idea more quantitatively. Fig. [7] displays 
the rms deviations between the ME parameters inferred from 
each line and the corresponding reference. The rms values 
have been calculated by dividing the x-axis in bins of size 
50 G, 9°, and 300 ms"^ The black and red lines refer to Fei 
525.02 and 525.06 nm, respectively. The solid lines indicate 
the results obtained from the spatially and spectrally degraded 
observations sampled at the five wavelength positions of IMaX 
(the corresponding maps are the ones in the right column of 
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Fig. 7. Rms diff'erences between the inferred ME parameters and 
the reference. The black and red lines stand for Fei 525.02 nm 
and 525.06 nm, respectively. The dotted lines correspond to the 
inversion of the degraded profiles but with full spectral sampling. 
The solid lines show the results for the degraded data sampled 
at the five wavelengths of the IMaX vector spectropolarimetric 
mode. 



Fig. [6]). The rms values refer only to pixels in which the Stokes 
Q, U or V signals of Fei 525.06 nm exceed three times the 
noise level, except for the LOS velocity where we included all 
the pixels. Clearly, the 525.02 nm line yields smaller deviations; 
typical values are 35 G for the field strength, 10° and 35° for 
the magnetic inclination and azimuth, and 50 m s"^for the LOS 
velocity. It is worth noting that the deviations caused by the 
limited spectral sampling are at least twice as large as those 
associated with instrumental degradation eff'ects (dotted lines in 
Fig. [71 see also the middel panels of Fig. [6]). 

In summary, the results show that both lines are suitable 
for magnetic field and velocity determinations based on ME 
inversions, although the deviations may be large for the field 
inclination and azimuth due to the weakness of the linear 
polarization signals. Among the two, however, the line at 
525.02 nm is to be preferred because of its larger signals and 
smaller deviations. 

3.2. Uncertainties of the IMaX spectral line 

The total uncertainties expected for the ME parameters inferred 
from the IMaX vector spectropolarimetric mode with the Fei 
line at 525.02 nm can be calculated directly by comparing the 
inversion results with the MHD models. Figure [8] shows the 
standard deviation of the diff'erences between the ME parameters 
and the corresponding MHD values at an optical deptlfl of 
logr = -1.5. Since an IMaX pixel is equivalent to four pixels 
in the MHD simulations, the latter have been binned by a factor 
of two in both directions to allow comparisons. For the magnetic 
field strength, the rms values are never larger than 100 G. For the 
field inclination, the uncertainty increases as the field weakens, 
reaching a maximum of ~ 30°. The rms values for the field 
azimuth are larger, on the order of 30-40°. In the case of the 
LOS velocity, the rms decreases for weak fields from 0.6 down 
to 0.6 kms"^ 



^ We have chosen the optical depth leading to the smaller rms 
differences. 
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Fig. 6. ME inversions of the Fei 525.02 nm line synthesized with the help of MHD models. From top to bottom: magnetic field 
strength, inclination and azimuth, and LOS velocity. The left panels represent the results from the data degraded by telescope 
diffraction and CCD pixelation. The middle panels show the maps retrieved from the degraded profiles after convolution with a 
6 pm wide filter and addition of noise. The right panels display the results obtained from the spatially and spectrally degraded 
observations with noise, sampled at five wavelengths points. Except for the reference, SNR = 1000 in all maps. 



There are three sources of error contributing to the values 
shown in Fig. (8) the uncertainties caused by the instrumental 
trade-off's (spectral smearing, wavelength sampling, and noise) 
as given in Fig. [71 the uncertainties associated with the spatial 
degradation, and the limitations of the ME approximation. The 
rms diff'erences displayed in Fig. [8] are largely dominated by the 
limitations of the ME analysis. 

A detailed study of th e uncertainties associated \y ith ME 
inversions can be found in lOrozco Suarez et a D (120101) . These 
authors analyzed noise-free synthetic observations computed 
from the same MHD models we use here, assuming very high 
spectral resolution and supercritical wavelength sampling^. They 
concluded that ME inversions provide atmospheric parameters 
that are averages of the actual parameters along the line of 



^ The study was carried out for the pair of Fe i lines at 630 nm. 



sight. Their results show that the uncertainties due to the ME 
approximation are larger than those associated with photon 
noise. 



4. Summary and conclusions 

We have presented an analysis of the reliability of Milne- 
Eddington inversions applied to high resolution magnetograph 
observations such as those delivered by IMaX. The analysis was 
indeed meant to help during the design and development phases 
of the instrument. In particular, the results of this work enabled 
the IMaX team to select the Fe i line at 525.02 nm instead of Fe i 
525.06 nm. Both lines could be used through the IMaX pre-filter, 
but the former permits more accurate inferences of the vector 
magnetic field. Fe 1 525.02 is also advantageous for the detection 
of the weakest linear polarization signals of the quiet Sun. 
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Fig. 8. Rms differences between the ME parameters inferred 
from the inversion of the Fei line at 525.02 nm as observed 
by IMaX and the corresponding MHD parameters at an optical 
depth of logr = -1.5. The different line shapes stand for 
different ME parameters: solid - magnetic field strength, dashed 
- field inclination, dotted - field azimuth, and dot-dashed - LOS 
velocity. 

Our analysis relies on synthetic Stokes spectra emerging 
from realistic MHD simulations of the quiet Sun. The spectra 
were degraded by taking into account all the SuNRiSE/IMaX 
instrumental effects, including telescope diffraction and CCD 
pixelation, smearing by the spectral PSF of the Fabry- 
Perot etalon, and coarse spectral sampling (the IMaX vector 
spectropolarimetric mode scans the polarization profiles at just 
four wavelengths within the line and one continuum point). 
Finally, the simulated IMaX observations were inverted in terms 
of Milne-Eddington atmospheres. 

Our main conclusion is that ME inversions can provide 
reasonably accurate atmospheric parameters (especially the 
magnetic field strength and the LOS velocity) from the 
resolved but coarsely sampled Stokes spectra acquired by 
SuNRiSE/IMa X. This conf irms the conclusions of previous 
works (e.g., iGraham et al.l 120021) . As a by-product, we have 
demonstrated the important effects of pure telescope diffraction 
on spectral profiles coming from the highly inhomogeneous, 
quiet-Sun photosphere. 
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